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Dynamical correlations in asymmetric infinite nuclear matter are investigates in a transport the- 
oretical approach. Self-energies due to short range correlations and their influence on the nucleon 
spectral functions are described in an approach accounting for a realistic treatment of mean-field 
dynamics and a self-consistently derived quasi-particle interaction. Landau-Migdal theory is used to 
derived the short range interaction from a phenomenological Skyrme energy density functional. The 
spectral functions in asymmetric nuclear matter are found to follow in their gross features closely the 
patterns observed previously in symmetric nuclear matter. An interesting sensitivity of dynamical 
self-energies and spectral functions on the momentum structure of the underlying interactions is 
found. 
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I. INTRODUCTION 
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£S| I Correlations beyond those described by the static mean-field are of genuine interest for the understanding of the 
dynamics in a nuclear many-body system. Although experience seems to confirm the prevalence of mean-field dynamics 
and the directly related quasi-particle picture a closer inspection reveals that nucleons in matter obey more involved 
rules. For symmetric nuclear matter this has been discussed extensively in the literature, e.g. in a recent review 
and also was considered in our previous work However investigations for asymmetric nuclear matter are rare. 

Some of the few cases are and 0. In symmetric matter and stable finite nuclei one finds typically a suppression 
of the quasi-particle strength by about 10%. Pure mean-field dynamics predicts that in nuclear matter the spin- and 
momentum states are uniformly occupied inside the Fermi sphere up to the sharp surface at the Fermi momentum 
. hp. Hence, the momentum distribution n q (k) for protons and neutrons (q=p,n), respectively, up to a normalization 
constant is described appropriately by a step function, n q (k) ~ 9(k F — k) where k F = k F for symmetric nuclear matter. 
As a matter of fact, in a more realistic description, going beyond the mean-field, this pattern changes. The nucleons 
cease to be in sharp energy-momentum states but acquire spectral functions of a finite width, thus resolving the strict 
on-shell relation between energy and momentum. Hence, the mean-field configurations only exist with a finite life 
time, determined by the strength of the coupling to 2h — lp and \p — 2h states, respectively. Correspondingly, in finite 
nuclei one observes a reduction of single particle spectroscopic factors when going away from the Fermi level. More 
dramatic effects are found in exotic nuclei. Leaving the valley of stability, the experiments H || and theory [ToL ITl| 
agree in observing strongly suppressed single particle probabilities already close to the Fermi level of neutron-rich 
nuclei. 

Despite the fact that strongly asymmetric exotic nuclei are of central interest for modern nuclear structure physics, 
the contributions of short range correlations to asymmetric nuclear matter and finite nuclei is still a rarely touched 
problem in nuclear physics. Short range correlations play an important role for the understanding of nuclear matter 
as for example seen in the observed high energy tails of nuclear momentum distributions and form factors. Also, short 
range correlations play a decisive role in understanding the structure of finite nuclei, particle production in heavy 
ion collisions and astrophysical objects like supernovae. Since supernovae are highly asymmetric in isospin we need a 
more detailed knowledge of correlations in asymmetric matter beyond ordinary mean-field interactions. 

Here, our main goal is to explore the interplay of mean-field dynamics and correlations in asymmetric nuclear matter. 
It is worthwhile to recall that in asymmetric matter isovector interactions are giving rise to strongly asymmetric Fermi 
energies for proton and neutrons. Hence, as an immediate consequence the available density of states of excited states 
changes, becoming increasingly different for proton and neutron configurations, respectively, with changing proton-to- 
neutron ratio. This a new aspect, clearly not accessible in the conventional approaches concentrating on symmetric 
nuclear matter. 

From our previous work we know that the spectral functions in symmetric nuclear matter are determined by 
interactions of extremely short range which are very well approximated by point couplings, reproducing almost 
perfectly results from rather involved many-body calculations, e.g. in 0] , both for cold and hot Q symmetric 
matter at all densities. It is therefore tempting to extend the same transport theoretical approach to asymmetric 
nuclear matter. The fully microscopic approaches to correlations in asymmetric matter use either only BHF self- 
energies from particle-particle ladders in the T-matrix approach |]| or include in addition hole-hole scattering 0, 
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, corresponding to 2hlp self-energies but neglecting the 2plh contributions for particle states. We go beyond 
the previous work and consider in this paper both types of dynamical self-energies, thus obtaining a consistent 
description of particle and hole states. However, this is achieved to the expense on treating the mean-field (BHF) 
parts phenomcnologically. 

A fully self-consistent description as indicated e.g. for symmetric matter in j|j is beyond the intention of the present 
paper. Therefore, long-range mean-field effects are included by a Skyrme energy functional by which we account for the 
above mentioned change in the Fermi energies of protons and neutrons, respectively. The basic transport theoretical 
relations are summarized in section [Tj] and are extended to the general case of asymmetric nuclear matter. We 
also discuss the influence of charge asymmetry on the mean-field, affecting both the effective mass and the effective 
binding potentials in asymmetric matter. Landau-Migdal theory is used to estimate the matrix elements for the 
ppp , nnn~ l , pnn~ l and npp~ l collisions by which the dissipation of spectral single particle strength is determined. 
Results of transport calculations for the spectral functions, self-energies and momentum distributions of protons and 
neutrons in asymmetric nuclear matter at various degrees of neutron excess are presented in section ITTT1 Finally we 
will close in section 17x71 with a summary and an outlook. 



II. SPECTRAL FUNCTIONS IN ASYMMETRIC NUCLEAR MATTER 

The underlying model for our calculations was presented in 0,01 and used for calculations in cold (T = 0) symmetric 
nuclear matter at saturation point density po and for finte temperature and high densities in At this point we 
restrict ourself to a short summary of the formalism with a few changes in respect to asymmetric nuclear matter. For 
a more detailed calculation see |3, Is L| • 

The two-particle-one-hole (2plh) and one-particle-two-hole (lp2h) polarization self-energies S > and S < for nucleons 
in nuclear matter are given by |13| : 

v-* f d 3 k 2 duj 2 d 3 k 3 dw 3 d 3 k 4 dcu 4 44 —- — ^ 
S ? = 9 }^ J -(2^)3 (2^)4 {2^^ $ (k + k 2 -h- fc 4 )|AV| 2 

xg>>{u 2l k 2 )gf(uj 3l k 3 )gf,(uj 4:1 k i ) (1) 

where g — 2 is the spin degeneracy. Compared to symmetric nuclear matter neutrons and protons occupy different 
Fermi spheres, which leads to different self-energies. In the asymmetric case one has to account for collisions with 
particles of the same kind (pp and nn) and for collisions with particles of different kind (pn and np). This explains 
the summation over the isospin q' in equation 

The one particle correlation function g q in equation Q can expressed in terms of spectral functions 

9q(u,k) = ia q {uj,k)9(uj q - u)), 

g>{cj,k) = -ia q (u,k)(l-0(uj q -oj)), (2) 

where uj q is the Fermi energy of the nucleons and 9{uj q — to) the step function to account the Fermi distribution in 
case of cold nuclear matter. The nonrelativistic spectral function for nucleons is found explicitly as, 



(<" - - - fc)) 2 + ir 2 ( w , k) : 



*( W > k ) = ~ S3 — v m/ ^ ^zr, rr^— v^r, ~ > ( 3 ) 



where Re£JJ e *(w, k) is the real part of the retarded self-energy, which can be calculated dispersively: 

Re^,k) = P [ d ^ T -4^± (4) 



The width T q (uj, k) is given by the imaginary part of the retarded self-energy, 



r,(w,fc) = 21ml]-* = i(E>(w,fc) - £<(w,fc)). (5) 



We introduce the asymmetry coefficient 



e = ^ (6) 
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which indicates the fraction of protons at a given total isoscalar nucleon density p — p n + p p . In a finite nucleus with 
charge number Z we have £ = Z/A. 

To see the influence of asymmetry we split the mean-field in proton and neutron contributions: 

= + T# = piu pp + P (l - i)u pn (7) 
= + E™/ = p(l - Z)u nn + pHu np . (8) 

E™^ and EJ^/ are generated by the proton background, which vanishes for pure neutron matter. E™ 1 / and ££Jj/ are 
generated by the neutron background. On the right hand side u pp and u pn are effective interactions averaged over 
the proton Fermi sphere for u pp and the neutron Fermi sphere for u pn in phase space. The same for u np and u nn . 

One can easily see in equation JSJ that there is no difference for neutrons and protons in case of symmetric nuclear 
matter (£ = 0.5) and one has not to distinguish between protons and neutrons. Going over to the asymmetric case 
the mean-field contributions split up, now protons and neutrons are distinguishable by their different dynamical laws. 
Absorbing the momentum dependent part of the mean-field up to second order in k into the kinetic energy term 

<*,,)■,**. + u;ff M ( 9 ) 

leads to an effective mass m* and an momentum-independent effective potential . Starting with symmetric matter, 
where the effective masses and effective potentials are the same, for p and n, respectively, the masses and potentials 
for nucleons of different isospin start to split up for asymmetric matter. As seen in equation the mean-field 
contribution plays an important role for the spectral function's pole structure. Hence, we expect a separation of the 
proton and neutron quasiparticle peaks. 

We account for mean-field effects in our model by means of a Skyrme energy density functional with parameters 
taken from the recent set of ref. |17| . By executing the variation with respect to the kinetic energy density r q and 
the nucleon density p qi where q = n,p one obtains for the effective mass and for the effective potential: 

^ = l + ^(^[ii(2 + .Ti)+t 2 (2 + a;2 )] /9 (10) 

+ ^[h(2x 2 + l)-t 1 (2x 1 + l)] P q) 

Ufl = i* [2(2 + x Q )p- 2(2z + 
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t 3 p a [2{2 + x 3 )p- 2(2x 3 + l) Pq ] 
at 3 p ff - X [(2 + x 3 )p 2 - (2x 3 + l)(p 2 p + pi) 



[t x (2 + xl) +t 2 (2 + x 2 )]T 

-[t2(2z2 + l)-ti(2xi + l)K (11) 
o 

where p is the total and r total kinetic density. The Skyrme parameters Xq, x\, x 2 , x 3l to, ti, t 2 , t 3 and a are given 
m tabled 

In the numerical calculations we use an iterative approach. Starting with an initial choice for the widths, i.e the 
imaginary part of the self-energy, we calculate the spectral functions, which then serve as input for re-calculating 
the self-energies. These results are used as input for the next step of an iterative calculation, thus improving the 
results. The calculation is stopped, when the spectral functions are converged to a given accuracy. In the language 
of Feynman diagrams this approach corresponds to the summation of the sunset diagram to all orders Q . 



A. The Short Range Interaction 



The former investigations were restricted to symmetric nuclear matter at saturation density. That highly symmetric 
situation simplifies the calculations because protons and neutrons are indistinguishable. Irrespective of the charge 
state, spectral functions and other observables can be related to the common Fermi energy of protons and neutrons. 
As pointed out before, in asymmetric matter the isovector mean-field changes the situation completely by introducing 
a relative shift of the proton and neutron Fermi surfaces, in addition to the changes from the differences in particle 
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TABLE I: The parameters of the Skyrme SLy230a [T3 interaction. The spin-orbit parameter Wo is included for completeness 
only. 



number. Also, by following the shift of the equilibrium density with increasing asymmetry, here we have to have 
information on the value of |M gg /| 2 at changing total number densities. However, from our former work we expect 
that only the short range parts of the full in- medium NN interaction will be relevant. 

An appropriate approach to investigate this aspect and, as an important side aspect, to test the consistency with 
the mean-field part of the model calculations we use Landau-Migdal theory [T^J . Given an energy density functional 
£{p) Landau-Migdal theory provides a way to calculate the corresponding residual 2-quasiparticle (2QP) interaction 
by second variation of £{p) with respect to the various spin and isospin densities |14l Hfij. Hence, for an intrinsically 
density dependent interaction the corresponding Landau-Migdal interaction includes various rearrangement terms 
which ensure to fulfill a number of important consistency relations |l4j . These relations would be violated if the bare 
G- matrix interaction would be used instead, as e.g. in 

For a Skyrme energy density functional the resulting parameters fo{p), fo(p), 9o(p) and g'o(p), characterizing the 
interaction strength and density dependence in the spin-isospin S = 0,T = 0, S = 0, T = 1, S = 1,T = 0, S = 1,T = 1 
s-wavep/i interaction channels, are in fact derivable in closed form, e.g. [l5j| . According to the structure of the Skyrme 
energy density functional the Landau-Migdal parameters are constrained to monopole and dipole components. 

However, these Skyrme values include short and long range interactions. The origin and nature of the various pieces 
cannot be identified directly. Here we are interested mainly in the short range parts. A simple but meaningful way to 
extract the short range parts is to identify the long range components with pion exchange. Consequently, we define 
the short range part of the interaction by subtracting from the interaction derived from the Skyrme energy density 
functional the (central part of) the pion exchange NN interaction. Hence, we use the central interaction obtained 
from pion exchange and include anti-symmetrization explicitly by means of the spin and isospin exchange operators 
P(j . t j respectively, thus leading to 

V^{k 1 ,k2)=-f^D v {k 1 ,k 2 )ffl-^Tl-T2{l-P cr Pr) , (12) 

which evidently can be arranged into 

V 7T {k 1 ,k 2 ) = V${k 1 ,k 2 )(a 1 -<J 2 ) S {T 1 -T 2 ) T (13) 

S,T=0,1 

Above, the momentum space pion propagator is denoted by 

DAh,k 2 )= * F{p 2 ) (14) 
(p z + m^) 

including a monopole form factor with a cutoff A = 800 MeV/c. The 3- momentum transfer is in the t-channel 

p = k 2 — k\ while in the u-channel we have p = Q = k 2 + k\. Numerically, we use / w ~ 0.075 which is the standard 
value for the pseudo- vector nNN coupling constant. 

Hence, we cast the energy density functional into the form 

E(p) = E s (p) + E„(p) (15) 

where E s = E — is the short-range contribution. The long range pionic part E % is given by a sum over the various 
spin and isospin transfer contributions as defined in ea !13l Formally E„ can be written as a sum over all spin-isospin 
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channels 



Mp)= \ E E /0/0 ^w-w^^)-fc 2 ) 

~,q>=p,n S,T=0,1 J 1 ' J 1 ' 



q,q'=p,n S,T =0,1 ' 



although in spin-saturated nuclear matter only the 5 = are non-vanishing. In symmetric, i.e. isospin saturated, 
nuclear matter also the T — 1 components will not contribute. The brackets indicate traces over spin and expectation 
values on isospin and the step functions constrain the momentum integrals to the proton and neutron Fermi spheres 
with Fermi momenta kp(q), q = p,n, respectively. 

Performing the variation with respect to the various spin and isospin densities we find in symmetric nuclear matter 
the set of Landau-Migdal parameters which in standard notation are 

/ W (p) = -4^AUMjrQo(l + ^) + / W (M (17) 

F F 

C\p) = +^N (2k F )^Qo(l + ^) (18) 

F F 



g£\p) = +4^ (fc F )J r Q (i + _g.) (,,), 



F 

1 I „ ... ml 



.'«(p) = 47^0^)^(^-8^00(1+2^)) (20) 



where Qo{z) = \ ln(f^f) is the Legendre polynominal of second kind for I = 0. A more detailed calculation is found 
in the appendix. In the S = 0, T = interaction channel an additional rearrangement term appears 

/M(fcir) = N (k F )V M (0,0)h(k f /m n ) (21) 

where (q, k) is the monopole component of the S — 0, T = ph-u-channel pion exchange interaction and 

, s (l + 2z 2 )ln(l + 4z 2 ) -Az 2 
^) = A ^3 1 ■ (22) 

The Landau-Migdal parameters are pure numbers obtained by normalizing the interactions to the level density at 
the Fermi surface. We use No(kp) — 2kpM* /(ttH) 2 , leading to 1/Nq = 154MeVfm 3 at the saturation point with 
kp = 1.33/m . 

The short range components are defined by subtracting the pion contributions, e.g. in the (S — 0, T = 0) interaction 
channel: 

/o (s) (p) = /o(p)-/o (7r) (p) (23) 

and / , g and g are defined accordingly. The £ = monopole Landau-Migdal Parameters for the full Lyon-4 
interaction and the short range parameters after subtraction of the pion parts are displayed in Fig|2 While the pion 
contributions are small in the (S = 0, T = 1) and the (S = 1, T = 0) channels, the (S = 0, T = 0) and {S = 1, T = 1) 
gain additional strength, being attractive in the first and repulsive in the second case. 

The short range components are obtained by subtraction of the pion part from the full Landau-Migdal parameters, 
e.g. in the S = 0,T = channel: /J(p) = fo{p) ~ fo(p) an d correspondingly for the other channels. The matrix 
element relevant for our calculations is defined from the pp, nn, np and pn spin-scalar and spin-vector amplitudes, 
respectively, 

UM = Mp) + (-) q '- q fo(p); (24) 

9*M = 9o(p) + (-) q '- q 9o(p) (25) 

by averaging over the spin degree of freedom. For the isospin projection quantum numbers we use q = ±g for neutrons 
and protons. Because of isospin symmetry we have f pp — f nn and f pn = f np , respectively. Hence, we define 

M q < q (p) = \(U q (p) + 3g q/q (p)) (26) 
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and averaging out the charge state dependence in a second step, 



M{p) = s i j-{\M pp {p)\ 2 + \M P n{p)\ 2 ) ■ (27) 

Results for M pp (p), M pp (p) and M(p) are shown in FigEI While the short range interaction among equal 
particles varies only slowly with density, the proton-neutron interaction shows a strong density gradient, remi- 
niscent to the known relations. At saturation density of symmetric nuclear matter, p = po = 0.16/m -3 we 
find |Af(po)| = 341Mey/m 3 . This value compares extremely well to the independently derived matrix element 
\M\ = 320MeVfm 3 determined in our previous investigations [H-EStH by adjusting the transport theoretical spectral 
functions to those from the many-body calculations of Bcnhar et al. [12J. This very gratifying agreement confirms 
the validity of our approach also for asymmetric matter. Moreover, with the derivation of the relevant interaction 
matrix element in a theoretically well justified approach we have gained predictive power in regions of total density 
p and for asymmetries £ = p p /p = Z/A < i which were inaccessible before. 

III. CORRELATIONS IN ASYMMETRIC NUCLEAR MATTER 
A. Interactions and Dynamical Self-Energies 

Spectral functions and the imaginary part of the self-energies have been calculated for nucleons in neutron-rich 
infinite nuclear matter in the transport theoretical approach outlined in the previous sections. Results will be displayed 
at the saturation density po = 0.16/m -3 of infinite symmetric nuclear matter although with decreasing values of the 
asymmetry £ the equilibrium point moves to smaller densities and finally disappears below a critical value of £. Our 
choice of the Lyon-4 parameter set [l7T | for the Skyrme energy density functional leads to an higher effective mass for 
protons than for neutrons in neutron-rich matter, i.e. m p (p,£) > m*(/9, £) for £ < i. 

The averaged scattering amplitude was chosen according to ea !27l thus accounting for the variation of the average 
matrix element with density. In all calculations, the real part of the self-energy, which plays an important role for 
analyticity , was included by means of using eq. (@J . By numerical reasons the spectral functions were calculated on 
energy and momentum grids (uj, k) with —0.5 GeV< uj < +0.5 GeV and < k < 1.5 GeV/c which obviously leads to 
a lack of information on the behavior of the spectral functions at energies outside the region covered by the numerical 
grid. For being able to calculate the ReS g by eq.(QJ, the imaginary parts were extrapolated into the regions of large 
energies by assuming Gaussian tails. 

The spectral properties of nucleons in asymmetric matter are best explored by keeping the total number density 
fixed to p — p eq . Then, the average matrix element is (\M qq '\ 2 )^ — 350 MeVfm 3 , being in agreement both with the 
derivation sketched above and reproducing the spectral functions of ref. |T^ . 

In the imaginary part of the self-energies dynamical correlations are reflected most clearly and directly. Results for 
ImE g = TjTg are displayed in fig^for protons and neutrons (q=p,n), respectively, at an asymmetry £ = j. The overall 
dependence of ImE g on energy is constrained by the requirement that the ground state is stationary. Theoretically, 
this is taken care off by using retarded Green functions, leading to Im£ 9 (wF(g)) = as seen in fig[S]and guaranteeing 
that the ground state is indeed stable as indicated by the vanishing width. 

In Fig[5] also results from a recent Brueckner Hartree-Fock (BHF) study of the Tubingen group Q are included. 
Both approaches agree qualitatively in the global energy dependence given by a strong increase of ImS g at energies 
away from the Fermi surface. But in detail, differences are seen to emerge on the quantitative level and especially in 
the hole sector: While our approach predicts long tails of ImS g extending to large values of \u — ojf\ and eventually 
levelling off by slowly approaching asymptotic values, the BHF calculations show a much more pronounced fall off 
at energies deep in the Fermi sea. Such a cutoff behavior is somewhat unexpected at first sight because the density 
of intermediate (2hlp) states still increases with increasing energy which should overcompensate the suppression 
introduced by the momentum structure of the interactions. The deviations in the particle sector (uj > ujp) are 
reflecting the differences in the interactions. The BHF calculations are based on the CD-Bonn NN-potential which 
seems to provide a larger amount of hard scattering at large off-shell momenta than our simplified description using 
a zero-range interaction, obviously not distinguishing between on- and off-shell processes. However, the results of 
ref.0 show that similar differences are observed in the BHF calculations when using different NN-potentials. This 
points to an interesting - and rare - sensitivity of the dynamical self-energies on the off-shell momentum structure of 
interactions. 

The striking similarity of our results to the full BHF calculations is a strong argument for the universality of the 
underlying dynamical effects. At first sight it might be surprising that the close relationship is obtained despite the 
fact that our approach does not include explicitly the tensor interaction. Tensor interactions play a crucial role in 
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BHF and also relativistic Dirac-BHF (DBHF) calculations 6]. The tensor interaction is of rank-2 in the spin and 
the orbital variables and as such has a vanishing expectation value for the ground state of spin saturated nuclear 
matter which we consider here. However, the tensor interaction contributes from second order on to the equation of 
state of symmetric and asymmetric nuclear matter. In the BHF and also the DBHF calculations the tensor parts 
introduce an especially strong proton-neutron correlation in the particle-particle channel 0, IS • Obviously, the 
on-shell part of that effect is accounted for by the phenomenological Skyrme functional but we cannot exclude that 
off-shell contributions might be missing. 

A slightly different situation is encountered in the dynamical self-energies. As seen from eq^ these self-energies 
are of at least second order in the residual interaction leading from the ground state to the intermediate excited 
states and back to the ground state. In these processes, the tensor interaction contributes in principle through the 
excitation of spin and spin-isospin modes in the particle- hole channel. However, since the corresponding particle- hole 
interactions are repulsive the spectral strengths of these modes are shifted to higher energies. Also, as seen in Fig|3 
the interactions in the spin-isospin channels arc relatively weak, inhibiting collective enhancements as observed in the 
spin-independent channels. Hence, the spectral properties close to the Fermi-surface are only weakly affected although 
contributions at higher energies and in the far off-shell regions cannot be excluded. Still, the close agreement with 
the BHF results of ref.jjj indicates that such effects cannot play a major role. 

In free space the tensor interaction plays a crucial role in the proton-neutron system as reflected by the binding of 
the deuteron. In the present context one could expect a correlation in the particle-particle channel. However, explicit 
calculations by other groups have shown that such correlations die out quickly with increasing density. After all, 
the bulk properties of nuclei are well understood by nucleonic degrees of freedom without the need to introduce a 
significant amount neither of condensed deuteron matter nor pn pairing in the isospin / — channel, to the best of our 
knowledge. The reason is that the in-medium tensor interactions and, accordingly, the deuteron-like particle-particle 
channels are found to be suppressed with increasing densities [2, IS- For example, the deuteron as a bound state 
disappears beyond the so-called Mott density |2C1 l2ll | . But as a reminder a pronounced pn correlation might survive 
at higher densities. As a side remark we note that the situation is different when investigating processes especially 
magnifying the in-medium pair channels like photo-absorption on nuclei above the giant dipole resonance where pair 
processes play a significant role (see e.g. |2^)- 



B. Proton and Neutron Spectral Functions 

In Fig. Elproton and neutron spectral functions at p — p eq are displayed for two different momenta and asymmetries 
£ = 4, i. In all cases, the quasiparticle peaks are clearly identified, as one would expect from symmetric nuclear 
matter. The main difference to symmetric nuclear matter are in fact the changes in the mean- field due to the 
isovector self-energy which translate into corresponding changes in the chemical potentials. With increasing neutron 
excess the protons are bound increasingly stronger while the effective neutron potential becomes shallower, leading 
to less binding. Thus, the number asymmetry induces via the isovector interactions a gap between the proton and 
neutron Fermi surfaces. Hence, both the shape and the location of the peak are influenced by the isovector interactions 
through their influence on the chemical potentials. 

Comparing the results for £ = j with those for symmetric matter, £ = i we find from Fig|B|that in asymmetric 
matter the proton hole strength functions are moved towards the Fermi level while the proton particle strength is 
moved away from the Fermi level. The neutron spectral distributions show the complementary behavior: the hole 
strength is shifted down into the Fermi sea while the particle strength is approaching the Fermi level. In neutron-rich 
matter the Fermi momentum for protons is smaller than in symmetric nuclear matter and for neutrons vice versa. We 
first consider the spectral functions in the hole sector, i.e. for states in the Fermi sphere k q < kF q . The hole spectral 
functions peak are close the mean-field on-shell energy position 



hk 2 ff h(k 2 -k\ n ) 

" = — ^U?f-u fq = —*SL . (28) 



'"" 2m* q " JH 2m q 



For a momentum of k = 31.3 MeV/c this leads to a proton quasiparticle peak at about uj v on — —28 MeV, while 
the neutron quasiparticle position is at = —78 MeV, as seen in FigHJ] Hence, the proton peak is clearly shifted 
towards the Fermi edge, while the neutron strength is repelled. Similar observations are made in the particle sector. 
For k = 312.5 MeV/c, corresponding to u>(p) — 35 MeV and ui(n) = 8 MeV respectively, the quasiparticle peak for 
the neutron is shifted towards the Fermi level, while the quasiparticle peak for the neutron is shifted away. 

Another observation from Figgis that the width of the quasi-particle peaks decreases when approaching the Fermi 
surface. This is in agreement with the general expectation that closer to the Fermi surface the life time of quasi- 
particle states increases strongly. The energy-momentum structure of r(w, k) is illustrated in Fig0|for two momentum 
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cuts in symmetric and asymmetric matter. From the behavior of T(u), fc) we conclude that in asymmetric matter the 
collision rates in the hole sector increase for neutrons while for protons they slightly decreases. For either type of 
nucleon the amount of correlations decrease for states above the Fermi level compared to symmetric matter. These 
results indicate that the real and imaginary parts of the polarization self-energies roughly scale with the density of 
states of the same kind of particles. Similar observation have been made in 0. 

The dependence of spectral functions in nuclear matter on density is investigated in FigHand FigEI Spectral func- 
tions for symmetric matter and at total densities between half and twice nuclear saturation density p eq = 0.16/m~ 3 
are displayed. The results are obtained by using the density dependent matrix element from equation 1271 In Fig[7| 
the spectral distributions of a nucleon well inside the Fermi sea are displayed. The evolution of the spectral functions 
with density reflect the changes in the chemical potentials. At p = \p eq a state with the chosen momentum is much 
closer to the Fermi surface than for p — p eq or p = 2p eq . This explains the variation in width and position of the 
quasi-particle peak seen in FigEI seen to be shifted into the Fermi sea with a considerable gain of width and high 
energy strength when the density increases. 

In FiglHl spectral functions for a particle state well above the Fermi surface are shown. Here, a complementary 
behavior is found. With increasing density the relative distance of the quasi-particle from the Fermi surface decreases 
which is reflected in a reduction of the energy shift and the width of the quasi-particle peak at higher density. In 
both cases, long spectral tails in the region below or above the Fermi surface, respectively, are found, indicating a 
considerable shift of single particle strength away from the quasi-particle peak. The overall features, however, are 
qualitatively similar to those observed previously in symmetric nuclear matter 0, Q although the results differ on the 
quantitative level. 



C. Proton and Neutron Momentum Distribution 



We finally consider the proton and neutron one-body ground state momentum distributions 

n q (k) = / duia q (uj 7 k) (29) 



where p q is the chemical potential for nucleons of type q = p,n. Results for symmetric and asymmetric nuclear 
matter are shown in Fig|§] There, the shrinking of the proton Fermi sphere in neutron-rich matter is clearly visible 
in terms of the diminished region of occupied proton momentum states. In asymmetric matter the structure of the 
momentum distributions follows otherwise closely the same pattern as found already in symmetric matter: Inside the 
Fermi sea the occupancy is reduced by about 10% and this fraction of the spectral strength is shifted into a high 
momentum tail with an almost exponential decline. At first sight it seem surprising to find that at high momenta the 
distributions approach the each other showing very similar magnitudes and slopes. This result is another illustration 
of the universality of the dynamical short range correlations which already was found in our previous investigations 
0,0, 01- The bulk parts of proton and neutron matter are seen to be contained in regions k < kp(q), where obviously, 
^f(p) < kp(n) in neutron-rich matter. The short range correlations produce a similar pattern as in symmetric matter 
by depopulating the Fermi sea by about 10% irrespective of the asymmetry. In fact in all cases the high momentum 
tails are almost identical, thus confirming the universality of the underlying dynamical processes. 

Integrating the proton and neutron momentum-distributions over a sphere with radius k in momentum-space and 
normalizing the result to the density we obtain the cumulative density index 

xJk) = — / dk'n q (k')k' 2 (30) 
PgTT Jo 

where x q (k) — > 1 for k — > oo. The quantity x q is particular well suited to display the fractional relative exhaustion 
of the given density up to a finite momentum k. For a system of non- interacting quasi-particles we expect x q (k) ~ 
(-^) 3 9(kp — fc) + 6{k — kp). In a system with correlations the Fermi surface will be visible as a discontinuity in the 
slope as seen in FigEI] The results displayed in Fig^Jshow very directly the significant difference in the correlations 
pattern of protons and neutrons: A large fraction of the proton strength is shifted into the high momentum region 
as indicated by the slow convergence of x p (k) towards unity. For symmetric nuclear matter we find from equation 1301 
that about 82% of the states in the Fermi sphere (fc < kp q ) are occupied, while in neutron-rich matter with £ = 4 
only 71% of the proton but 87% of the neutron states inside the Fermi sphere are occupied. This observation is on a 
qualitative level in agreement with BHF calculations Q . Interestingly, we have obtained the asymmetric depletion of 
proton and neutron occupation without using explicitly a tensor interaction. 
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IV. SUMMARY AND OUTLOOK 



Dynamical correlations in asymmetric infinite nuclear matter were investigated in a transport theoretical approach, 
thereby extending our previous work into a new regime. It is worthwhile to emphasize that only few studies on 
dynamical correlations in asymmetric exist. In asymmetric matter isovector effects are strongly enhanced with in- 
creasing asymmetry. On the mean-field level this was taken into account by using a modern Skyrme energy density 
functional from the Lyon group |17| with parameters adjusted both to finite nuclei and neutron stars. In this sense, 
we describe the evolution of mean-field dynamics with asymmetry and density in a realistic model. In particular, the 
saturation properties of infinite nuclear matter are well described. Using Landau-Migdal theory we could derive from 
the Skyrme functional the corresponding effective quasi-particle interaction in a self-consistent way. An approach was 
presented which allows to extract the short-range interactions by subtracting the pion exchange contributions. 

The overall features of self-energies and spectral functions in asymmetric nuclear matter resemble those found in our 
previous work on short range correlations in symmetric matter. In particular, the present investigation confirms our 
previous conjecture on the universal character of short range correlations in infinite nuclear matter. An interesting 
observation in the present context is the sensitivity of the spectral distributions on the momentum structure of 
self-energies and interactions. 
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APPENDIX A: LANDAU-MIGDAL PARAMETERS FOR THE PION 



Second variation of the energy density in equation (|16[) with respect to spin and isospin densities leads to 
d 2 E n (p) 

S.T=0.l 

" ' " < 0?i •<? 2 ) 6 '(ri -t 2 ) t >] 



— — -> = J2 [VsTikf^kfiq'M^-azfin-rzf (Al) 

OpqOPqi 



i^W rffcfc 2 e(fc/( g ")-fc)|- v ST 



where the brackets indicate ground state expectation values. From which we now identify the Landau-Migdal param- 
eters: 

/ W = N (k f )V oa + f r (k f ) (A2) 
/'« = N (k f )V 10 (A3) 

gW = N (k f )V 01 (A4) 
g 'M = N (kf)V n (A5) 

where f r (kf) is an rerrangement term, which appears for spin-saturated symetric nuclear matter only in the S = 0, 
T = channel. Expansion of the in-medium scattering amplitude in a Legendre series, which is common, leads to 



1 f 1 

./f )= 2/ Pl(cos6)r(cos0)dcos9. (6) 



The integration over the angle 9 is done using the Newman formula 

1 f 1 Pl(x) 



(7) 
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where the Qi(t) are Lendendre polynominals of the 2nd kind. Using equation Q together with our expressions for 
the Landau-Midgal parameters leads to 

t ] = -^0(^)^*1 + ^) + ff\k F ) (8) 

F ' F 

//W = 4^iVo(M^QKl + ||) 0) 

F F 

I \ 3 f 2 T77 2 

= 4^0(^)^(1 + ^-) (10) 

ft' (7r) = 9^0-4^(^)^,(1 + ^1). (11) 

where the / = components in the S = 1, T = 1 channel includes in addition the direct (i.e. Hartree-type) 
contributions 

g' H = 4^V (fc F )4" (12) 
For the additional rearrangement term in equation (jSJ) we get 

f {r \k F ) = iv (fef)i / ' dkk 2 e(k F -k)^- v^\k,k') (13) 

where V^f (fc, fc') is the monopol component of the S = 0, T = u-channel pion exchange interaction. 
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FIG. 1: Equation of state of infinite nuclear matter at various asymmetries £ = p p /p = Z/A. Results obtained with the Lyon-4 
parameter set are shown. 
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FIG. 2: Landau-Migdal parameters in symmetric nuclear matter for the full Lyon-4 interaction (upper panel) [rUl and the 
pion-subtracted short range interaction (lower panel). 
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FIG. 3: Dependence of the average short range interaction matrix elements \M qq ' and of the average matrix element M , ca l27l 
on density in infinite nuclear matter. 



15 




FIG. 4: The nucleon width for symmetric (dashed line) and asymmetric nuclear matter (full line). The upper graphs refer to 
protons the lower ones to neutrons. 
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FIG. 5: The imaginary part of the nucleon self-energy for momentum k = OAkf and asymmetry £ = 0.25 at a density p — 0.17 
l//m. The dotted line represents the BHF results from Q. 
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FIG. 6: The nucleon spectral function for symmetric and asymmetric nuclear matter. The dashed line refers to symmetric and 
the straight line to asymmetric nuclear matter with an asymmetry of £ = 0.25. The upper graphs show the results for protons, 
the lower ones the results for neutrons. 




FIG. 7: Spectral functions of a nucleon well inside the Fermi sea at density p = \p eqi p = peg and p = 2p, 
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FIG. 9: The momentum distribution for protons (upper graph) and neutrons (lower graph). The straight line is the result for 
symmetric nuclear matter the dashed line for asymmetric nuclear matter (£ = 0.25) 
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FIG. 10: The the cumulative density index for symmetric nuclear matter (full line) and for protons (dashed line) and neutrons 
(dotted line) in neutron-rich matter with £ = 0.25. 



